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Abstract 

Radiative transfer equation in plane parallel geometry and Kompaneets equa- 
tion is solved simultaneously to obtain theoretical spectrum of 1-125 KeV 
photon energy range. Diffuse radiation field is calculated using time- independent 
radiative transfer equation in plane parallel geometry, which is developed 
using discrete space theory (DST) of radiative transfer in a homogeneous 
medium for different optical depths. We assumed free-free emission and ab- 
sorption and emission due to electron gas to be operating in the medium. 

/ dn \ 

The three terms n, n 2 and ( — — J where n is photon phase density and 

Xk = (tt^~) > m Kompaneets equation and those due to free-free emission 
V kT e / 

are utilized to calculate the change in the photon phase density in a hot 
electron gas. Two types of incident radiation are considered: (1) isotropic 
radiation with the modified black body radiation I MB [1] and (2) anisotropic 
radiation which is angle dependent. The emergent radiation at r = and 
reflected radiation r = r max are calculated by using the diffuse radiation from 
the medium. The emergent and reflected radiation contain the free-free emis- 
sion and emission from the hot electron gas. Kompaneets equation gives the 
changes in photon phase densities in different types of media. Although the 
initial spectrum is angle dependent, the Kompaneets equation gives a spec- 
trum which is angle independent after several Compton scattering times. 
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1. Introduction 

Compton and inverse Compton scattering play an important role in the 
processes of emission of X-ray spectrum which has been observed in many 
compact astronomical objects. Comptonization of X-rays has been dealt 
mostly through Fokker-Planck approximation (i.e, electron temperatures and 
photon frequencies are not high compared to the rest energy of the elec- 
tron) [2] [3] [4] [5] and others. Primary non thermal X-ray and 7-ray emission 
may reprocessed and end up in different energy band [6]. Several X-ray ob- 
servations obtained by Ginga do not confirm the power law spectral shape [7] , 
[8]. Spectral hardening have been noticed from Ginga satellite observations 
[7] [9] [10] and [11] of Seyfert galaxies and Cygx-1. 

Analytical approximations of Compton reflection have been developed by 
[12], on the basis of the separation of spacial and energy transport in the 
transfer equation. This approach have the advantage of facility to explore 
the model parameters. Later, many authors [13] [14] [15] [16] and [17] to 
implement Monte Carlo methods to treat several aspects of Compton re- 
flection models by taking more detail of the physics and geometry of the 
problem. Burigana [18] studied the reflection problem with semi analytical 
transfer equation and obtained the solution of the transfer equation. Angle 
dependent Compton reflection of X-rays and 7-rays has been studied by [19]. 
Studies on Compton reflection of X-rays and 7-rays by [12] [20] [21] and oth- 
ers in shaping the emergent X-ray and 7-ray spectra. Czerny & Elvis [22], 
Zdziarski et al [23] and Haardt & Marasachi [24] studied the inverse Comp- 
ton scattering in producing the high energy photons in some of the compact 
sources in which low energy photons gain energy from hot electrons. The 
geometry and distribution of matter in these objects would make the process 
anisotropic [25] [6] [26]. Haardt & Marashi [24] and Haardt [27] studied this 
process in plane parallel geometry. Nishimura, Mitsadu & Itoh [28] studied 
the same problem but without angular dependence with isotropic scattering 
in the laboratory frame with kT e < 150 KeV. 

The escape probability and the energy distribution of the single scattered 
photons are analytically calculated using Mellin integral transformations. 
Miyamoto [29] estimated the photon escape time distribution and effect of 
Comptonization on the spectrum. Katz [30] obtained space distribution of 
the photons and spectrum changes by solving Fokker-Planck equation from 
Shapiro et al [31]. The transfer of radiation estimated by using Monte-Carlo 
method [32]. Another characteristic of X-ray spectra is its variability over 
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time. The time delay due to photon diffusion and the energy release and mul- 
tiple scattering of X-ray photons with thermal electrons due to comptoniza- 
tion would change the nature of the emergent X-ray spectra. If kT e > hv 
then the photon gain energy and if kT e < hv then the recoil effect will re- 
sult in decrease in photon energy. Kompaneets equation [33] allows us to 
find the time evolution of a given initial spectrum due to comptonization in 
homogeneous medium. Illarionov & Syunyaev [34] obtained few analytical 
solutions with a given initial spectrum in the soft X-ray region for continuum. 
Sunyaev & Titarchuk [25], Sunyaev &Titarchuk [35] calculated few solutions 
analytically for the time evolution of a given initial soft X-ray spectra with 
Kompaneets equation and polarization of hard radiation. The presence of 
electrons with random velocities would increase the energy of the photons 
Av V 2 kT 

by ps — — ps — |. This is the essential information that the Kompaneets 

v c l mc 2 

equation Kompaneets [33]. Wehrse & Hof [36] have made line calculations 
using inverse Compton scattering. 

The problem of the structure of a model atmosphere and the spectrum 
formation under the influence of very hard X-ray external irradiation is not 
satisfactorily solved in the present literature. Recent papers with more 
general descriptions of Compton scattering were published by many au- 
thors [37] [38] [39] [40] [41] [42] [43] and [44] on X-ray irradiated model stellar 
atmospheres. Majczyna et al [45] had presented a model which includes 
set of plane-parallel model atmosphere equations for a hot neutron star. 
Suleimanov et al [46] [4 7] examined the effects of Compton scattering on the 
emergent spectra of DA white dwarfs in soft X-ray range. They adopted two 
independent numerical approaches to the inclusion of Compton scattering 
in the computation of pure hydrogen atmosphere in hydrostatic equilibrium. 
The Kompaneets diffusion approximation formalism is used in one case, cross- 
sections and redistribution functions of [48] is used in another case. Chluba 
& Sunyaev [49] presented analytic solution of the Kompaneets equation for 
low frequencies. They noticed that multiple scattering of photons by hot 
electrons lead to brodening and shifting the spectral features. Exact station- 
ary solutions of the Kompneets kinetic equation is given by [50]. He also 
showed that the photon input and output points always corresponds to fi- 
nite frequencies. Procopio & Burigana [51] developed numerical code for the 
solution of Kopaneets equation and discussed accuracy and applicability. 

In this paper diffuse radiation field is calculated to estimate the reflected 
and the emergent radiation out of a plane parallel slab and also to calcu- 
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late the comptonization of the photons with energy between lKev-125Kev. 
Radiative transfer equation in plane parallel approximation is solved with 
free-free absorption and emission and emission from hot electronic gas in the 
medium. Photon phase density is estimated from the specific intensity de- 
rived from the equation of transfer for both cases isotropic and anisotropic 
incident radiation. Kompaneets equation is used to calculate time evolu- 
tion of photon phase density. At present we are not aiming to explain the 
observational results. 



2. Description of the Theory: 

Kompaneets equation: The photon energy in an electron gas with ran- 
dom velocities, would on the average increase by amounts approximately 

ps — Ps — %■ where k, T e , m, c, v are the Boltzmann's constant, tem- 

v C l mc 2 

perature of the electron gas, mass of the electron, velocity of light, frequency 
of the photon respectively. The photon phase changes over time due to scat- 
tering. This problem is well described by the Kompaneets equation which 
gives the time evolution of a given initial photon spectrum due to Comp- 
tonization and is given by [33] [1], and [25], 



dn 
dy 



1 d 

xl dx k 



n + n z + 



dn 
dx 



(1) 



where 



hv 



(2) 



n 



2hv 3 



(3) 



V = — e -a T n e cdt (4) 
Jo mc 

where I is the specific intensity, h is the Planck constant, a T is the Thom- 
son scattering cross section and n e electron density. The n 2 term in equa- 
tion (1) describes the induced Compton interactions, ie., the energy being 
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transferee! from hotter photons to cooler electrons as a result of which the 

photon energy is reduced [52] and [53]. The term is important in the low 

frequency range, the term containing n describes the drop in the photon en- 
Az/ hv 

ergy m in the photon energy through scattering which goes to the 

v mc l 

( dn \ 

heating of the electrons. The term containing — — represents the diffusion 

V ox k J 

of photons mainly to increase the photon energy with consequent cooling of 
electrons. 

Free-Free (Bremsstrahlung) processes influence the frequency redistribu- 
tion of the photons [33] and [34] and is given by 



where the coefficient a is given by 



mc 

and 



kT 

\a T n e c (6) 



k = 1.25 x lO- 12 g(x k )n 2 e T- 3 - 5 

= k g(x k ) (7) 

g(xk) being the Gaunt factor given by 

{1 for x k > 1 
V3 } /2.35x (8) 
In I for x k < 1 
7T V x k J 

Therefore the complete equation which includes both free-free and Compton 
processes is given by 

dn /dn\ /dn\ 
dy \dyJc \dy/ff 

ke~ Xk 

+ 



ax\ 



■{l -n(e x " - l)} (9) 
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Equation (9) is solved numerically on discrete mesh of frequency points. We 
obtain the spectrum of occupation numbers of photons in the phase space as 
n recursive relation given by 



where 



n i+2 = (p t A^j ni{ax - + n i+1 (l - a 2 ) 
-p t ^A 10 + M^rii+i)} 

on = Pt ( - 2Ai + A 4 + 4A 5 - A 6 + An) 



(10) 



(11) 



a 2 = p t (2Ax + A 4 + 4A 5 - 2A 6 - An) 



(12) 



M (ni, n m ) = n? (A x - A 4 ) + n^ +1 ^! + A 4 ) + 2A 1 r W+1 (13) 



where 



Ai = \(xl +1 + xl); A 2 = A\- ^=(Ax fc ) 1 ; 
A 4 = A 2 A 3 ; A 5 = AiA 3 ; A 6 = A 2 A§; 

A 7 = A?; A 8 = e- Al ; A 9 = e Al = (a 



A w = FA h A x 3 ; A u = ^A 10 (a 9 - l) 
ft = At.P; P = = 3.37 x 10~ 24 T R n P 



mc 



F = 3.7147 x 10 n #(x)n e T e 



-1-4.5 



_ 5.92648 x 10 



9 



hu 



(14) 
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2.1. Solution of Radiative Transfer Equation 

We obtain the initial spectrum of n in the equation (3) by solving the 
equation of transfer for I, the specific intensity is given by [54], [55] and [2] 
and others 

\dl v dl v . 1 (dl u \ 

c at or c\ at / scat 

where I v (r,l,t) is the specific intensity of radiation in 
ergs cm~ 2 S~ 1 H z ~ 1 str~ 1 , at position r, time t, and frequency v traveling 
in the direction given by the unit vector I. In the above equation, it is 
convenient to use the photon energy e = hv rather than frequency. (I u ) 

can be written as /(e) = -j- (KeV cm~ 2 S~ 1 KeV~ 1 str^ 1 ). The quantity 

na v = S ^^ n,jO,j{y) is the absorption coefficient, rij is the number of absorbers 

3 

of type j and o~j(v) is the total cross section of true absorption of photon 
energy. All effect of electron scattering are included in the term 

V ijt J scat 

These were studied by [54] which are valid for IKeV < e < lGeV for which 
the photon occupation number n is less than unity which is valid for many 
X-ray and 7-ray sources. If n « 1 then the induced scattering is neglected. 
In a cold electron gas for e > IKeV, [54] gives, 

(iL =cne / dw ®K r ' r ' t ) -K r ' r ' *)] (16) 

pc p 

correct to the lowest — and — where p is the momentum of the targeted 

e mc 

charged particle in the observers frame of reference, T = - is the photon 

c 

momentum before scattering and V is the momentum after scattering [55]. 
These are related as 

r 1 

r=r^ < 17 > 



1 



where 6 is the angle between the vector moments of T and T and — is 

ail 

the differential cross section for Compton scattering in the observer's frame 
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where the plasma is at rest. This is given by Klen-Nishana formula from [56] 
and [57] 

da i 2 /r'\ 2 /r' r 2 \ 

c 2 

where r$ = — - = 2.85 x 10 13 cm is the classical radius of the electron. 

mc 2 

For an incoming photon energies e' small compared to the mc 2 , which is 

/\f f / dn\ 

seen from the equation (17) that << 1. The collision term ( — ) 

e \ at J scat 

is then given by Fokker-Planck type expression derived from (16) [58] and 

[59]. The resulting equation of transfer is given by (to the lowest order of 

* = ^ = ^) [55], 
mc 2 mc 2 

ld6I x T dl x , w . d / 



c dt 



+ I— + n e a T cr(x)I x = j x + a T n e x— [xl x J (19) 



where n e orO[x) is mass absorption coefficient, j x is the emissivity 

in ergs cm~ 3 s^ 1 H z^ 1 str^ 1 . 
And 

a{x) = a s {x) + a a {x) + a[ f (x k ) (20) 
where for hv « -mc 2 [56], the quantities a x (x),a a (x) are given by 

a s (x k ) = 1 - 3x + ^x 2 - 28x 3 + ^x 4 - . . . (21) 



42 2 147 o 1616 4 

aJx) = x x 2 + x 3 x A + ... (22) 

v ; 10 10 35 V ; 

In the present analysis we restrict our study to the energy range of 1 KeV-125 
KeV (hv < -mc 2 ). The free- free absorption coefficient cr^J (x) is given by [1] 

a f J(x k ) = 1.79 x lO- 18 ^ 1 Z 2 n lX - 2 g ff (23) 

where z is the charge and rij is number density of ions and g~ff is the Gaunt 
factor which is of the order of unity. The quantity j(x) is the emissivity is 
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given by 



j(x) = -B&(x k ) + B e (x,P) 



+— / R(x,x')I(x',fj,')dx'dfj,' 
4vr J 

where B^(x) is the Bremsstrahlung emission given by [1] 
B£(» = 6.8 x 1(T 38 Z 2 ni n e e "& g ff 



(24) 



(25) 



B£(i/) = 6.8 x 1(T 38 Z 2 mn e x 



T e 2 exp 



5.926 x 10 



9 , 



-:r 



Further the quantity B(x,/3) is given by [1], [27], and [28] 



B e (x,0) 

n e a T ld(x - x ) 

4~f 2 f3 2 e 



x 



x 



^ * 1 ~ ^ / x ° / 1 

^ /or TT^ < - <1 

xq 1 + /3 



(1 + /?)--(!-/?) M 1<^< 
x x 1 

otherwise 



(26) 



(27) 

where — xo) is the monochromatic specific intensity, (3 = — and 7 = 

1 c 
(1 — [3 2 )~ ■ B e (x,/3) becomes effective only at higher temperatures as the 

energy is transformed from hot electron gas to the cold photon gas. We need 

to study the redistribution of energies at different frequencies. This is given 

by [55], and [56] 



R(x, x' 

3 OTTi e 

8~^~ 



x' X / 1 1 \ / 1 1 \ 

-+--2 -"-+--- 

X X \x' X/ \X X/ 



2 n 



(28) 
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where x, x' are the energies of photons before and after the scatterings. Equa- 
tion (19) is valid for hv « mc 2 [55]. The last term on the right hand side 
of equation (19) is the same as the term containing n in the Kompaneets 
equation (1). Therefore equation (19) can not fully represent the diffusion 
of photon energy which is adequately done by Kompaneets equation as this 

contains the — — and n 2 term also. We need to solve equation (1) and (19) 

simultaneously without the last term on the RHS of equation (19). Therefore 
the time independent plane parallel equation of radiative transfer is written 

as 

dl(x,n,z) 
u ; h n e a T <7{x, 1 jl ix, fi, z) 

dz 

= a a (x)B e (x, (3) + aff [ X ' Te) Bf/(x, T e ) 

poo n + 1 

+ -1 / / R(x,x')I(x',^W (29) 
»tt Jo J-i 

where // is the cosine of the angle made by the ray with normal to the plane 
parallel layers, important steps of the method is given Appendix and 



a" = n e a T a a ff (x, T e ) = n € a T a f J . (30) 

3. Result and Discussion 

3.1. Boundary conditions 

We shall give the boundary condition of the incident radiation as follows: 
Case 1: 

I( T = T,n) =1 mb (31) 

and 

I(t = 0,a*) = (32) 

where 

I^B = g_4 x \^T}p^ n x\c^ [c Xk - l) 2 (33) 
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is the modified black body radiation given in [1]. Here T e is the temperature, 
p is the density, cjff is the Gaunt factor given by 



1 for x k > 1 

Sjf) = { VS^x < , (34) 

7T \ Xh. / 



hv _ 5.92647 x 10 9 /«/ 

£fc — : — — X] X — - J 



kT e T e mc 



Case 2: 



I(r = T,/i) = I 



MB 



X 



6 

(3 + 26) 



(36) 



with 6=2 and 



l(r = 0,^) = (37) 

in this case we set B e (x,/3)=0 but free-free emission is included in the 
medium. 

The solution of transfer equation is developed using DST and obtained the 
solution. The stability and the accuracy of the method is checked through the 
following; We divide the medium into a number of " Cells" whose thickness is 
less than or equal to the critical (T crit ). The critical thickness is determined 
on the basis of physical characteristics of the medium. T crit ensures the 
stability and uniqueness of the solution and other steps are mentioned in 
Appendix. By considering above, twenty trapezoidal points for energy, four 
angle points and fifty plane parallel homogeneous layers are considered in 
our calculations. As there are many parameters to be studied we choose few 
representative of these to highlight the physical processes. We have chosen 
temperature T e = 5 x 10 9 K and 10 11 K, n e = 10 15 cm _3 ,r = 1 with the 
initial conditions given in equation (31) to (37). We assumed a purely ionized 
hydrogen gas so that p in equation (33) becomes n e m p where m p is the mass 

of the proton. The energy points(x = — ^r) are chosen from 0.0002 to 0.24 

mc 2 

on trapezoidal quadrature corresponding approximately to 1 KeV and 125 
KeV respectively. The four angle points are chosen on the Gauss-Legendre 
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quadrature on (0, 1) with ^=0.6943, // 2 =0.33001, /i 3 =0.66999, /j 4 =0.93057 
and the corresponding weights are Ci=0.1739, 6*2=0.3261, C3 = C2, C4 — C1. 
Case I: The incident radiation field represented by I MB the modified black 
body radiation (ergs cm -2 S -1 Str" 1 ), is plotted in figure 1 for energy points 
at temperature 5x 10 9 , 10 10 , 10 11 K to see how a photon of given energy 
would get transformed through the medium for different temperatures. 

In figure 2, We plotted the variation of I~ , I + and their corresponding 
photon phase densities n + and n~ with energy for the parameters are shown 
in the figures (see equation A. 20, A. 21, and 3) for /z 4 when initial intensity 
is given at I~(N = 50) = 1 and free-free emission is not included. Figure 
2(a) shows the emergent intensity and figure 2(b) gives reflected intensity. 
Figure 2(c) and 2(d) give their corresponding photon phase densities N=50, 
25, 1 corresponding to r =0.5 and r =0. We notice that, a small fraction of 
emergent intensity is reflected at N = 50 and emergent intensity I~ increases 
marginally at higher energy points. 

In figures 3(a, b c, d) we show the variation of I~(fj), I + (fJ<), n~(fj), 
n + (y) across the medium for the photon energy 100 KeV. at N = 50) 

increases as it acquires more free-free emission in the medium. I~(fii) shows 
maximum values at about N=40 and falls slowly towards N=l. / + (/ii) > 
I + {^2) > I + (fJ-3) > I + {^a) at all layers. The intensities fall sharply from 
N=10 reaching minimum at n — 1. 

Figures 4(a, b, c, d) contain the information of I + , n~ , n + with respect 
to e at T e =10 10 K. The incident intensity I~ at N=50 is given as I MB and 
free-free emission from the medium is assumed. I~(e) reduces toward higher 
energy points although there is a slight increase at e psI-10 KeV. I~ at N — 1 
that is at the emergent point is considerably reduced compared to the that 
at N — 50. It appears that free-free emission at T e =10 10 K is not effective in 
contributing to I at higher energies. The reflected intensities I + at N = 50 
is higher than that at N — 1 which shows that more energy is reflected at 
iV = 50 than that at N = 1 at the emergent side (r = 0). At = 50 , 
I is greater than I + by a factor of 10 4 while at N — 1, I~ is larger by a 
factor of 10 at lower energies and I ~ I + at higher energies. The n~ and 
n + given in 4(c) and 4(d) reflect the same variations shown figures 4(a) and 
4(b) respectively. 

Figures 5(a) and 5(b) contain the transfer of an 8 KeV photon energy 
through the medium from iV = 50 to I. I~(8KeV, /xi) falls rapidly between 
iV = 50 and N — 1, that is from r = r max = 1 and r = 0, changing by a factor 
of 10 4 while I~(8KeV : /i 2 ) change by a factor of 10 and I~(8KeV, /i 3 , yU 4 ) 
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change only by a factor of approximately 2. The situation with I + is different 
in that I + (8KeV, Hi) change is not as steep as those of I + (8 KeV, ^2,^3, ^4). 
All I + (8 KeV) fall rapidly in the extremely outer layers from N w 10 to 
N — 1. A similar variation can be seen in figures 5 (c) and 5(d) for n~ and 
n + corresponding to I~ and /+. 

Figure 6(a, b, c, d) shows the enhancement of photon phase space density 
given by the recursive relation given in equation (10) for temperature T e = 
10 10 i^, r = 1 and the emergent intensities I~(N = 1), with I~(N = 50) = 
I MB and free-free emission. We have calculated the enhancement of rC for 
different values of At(see equation 13). In the recursive relation of (10), 
we choose the first two energy points and n 2 from I~(N = l,e = ex) 
I~(N = l,e = e 2 ) respectively. The photon phase space densities in figure 
6, (for /ii and ^4)) rise sharply from At = 100. It is interesting to note that 
n _ (/ii) are much smaller than those of n~(fj,4) up to e 50 KeV. Above 50 
keV, n~(^4) and n~(/ii) coincide and the phase space densities become angle 
dependent. 

Case II: Calculations are done with anisotropic boundary conditions given in 
the equations (36) & (37) due to [25]. Figure 7(a, b) represents I~ and I + 
for r = 1 for the boundary condition 

and 

r( T = o,/i) = 

with 6=2 and T = 10 10 K. We assumed free-free emission in the medium. 
No emission from the hot electron gas is given in the medium. There appears 
to be a maximum at about 8-9 KeV and a steep fall in the energy range of 
e > 10 KeV, More energy seems to be emerging through ^ than through 
/i 4 , particularly at N=l. The reflected I + is considerably reduced and is flat 
irrespective of energy range. Figure 8 (a, b) give the result of Fs for r = 1 
with the initial condition 

I~(r = w,/i) = I MB x 2^ Jn(^) 
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e(KeV) 

Figure 1: I MB (in ergs cm~ 2 S^ 1 H Z^ 1 Str~ 1 ) versus photon energy e in KeV, for tem- 
peratures T e = 5 x 10 9 K, 10 10 K, 10 n K 

4. Conclusions 

Using discrete space theory of radiative transfer comptonization spectra 
is calculated in a plane parallel atmosphere with isotropic and anisotropic 
as a incident radiation. The accuracy of the solution is assumed by observ- 
ing 1) strict radiant flux conservation 2) non-negativety and continuity of the 
solution across (— 1 < fj, < 1) by taking the step size for the democratiza- 
tion. We estimated the photon phase densities from the specific intensities 
which is calculated from radiative transfer equation and the time evolution 
of photon phase density by employing Kompaneets equation. There is sub- 
stantial differences are noticed in the nature of emergent intensity variations 
with free-free emissions and absorption cases. We also studied when a pho- 
ton of given energy would get transformed through the medium and noticed 
that when transfer takes place radiation field becomes angle dependent. We 
noticed that the initial spectrum is angle dependent and the Kompaneets 
equation gives a angle independent spectrum after many Compton scatter- 
ings. We are plan to extend the above method for spherically symmetric 
approximation to study the stellar atmospheres in X-ray binaries and accre- 
tion discs atmospheres in active galatic nuclei(AGN) etc. 
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Figure 2: (a) I~ the emergent intensity (in the direction of t — > 0) at different layers at 
N=50 (r = T max ), N~25(t — 0.5) and N=l (r = 0) for angle m for temperature 5 x 10 9 
K. (b) I + the reflected intensities at N=50, 25, and 1. The initial incident intensity is 
I~(N = 50) = 1. (c) and (d): shows the photon phase densities n~ and n + corresponding 
to I~ and / + given in figures (a) and (b) respectively 
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10 u 



10 c 





Figure 3: The transfer of photon with e=100 KeV energy across the medium from N = 50 
to N = 1, in a medium with free- free emission and I~(N = 50) = I MB as the initial 
intensity. 
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Figure 4: Same as those given in figure 3, but with r e = 10 10 i^ 
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Figure 5: The transfer of photon with e=8 KeV energy across the medium from TV = 50 
to N = 1, in a medium with free- free emission and I~(N = 50) = I MB as the initial 
intensity for T e =lQ w K 
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Figure 6: The photon phase density n~ corresponding to I" is given versus e. The 
enhancement of n~ for 5t=l, 2, 5, 10, 100 is shown. The incident radiation is I~(N = 
50) = I MB in a medium with free-free emission 
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Appendix A. 

Equation (29) can be solved numerically by using discrete space theory 
of radiative transfer [60], [61], and [62]. We discritise the equation (29) on 
the [x(i) — z(n) — /x(m)] discrete mesh and this becomes 
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T represents the transpose and /i's and C's are the roots and weights of angle 
quadrature and 
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Similarly R H , R + and R . In case of static medium, 

R++ = R+- = R-+ = R (A.10) 

The subscript n + \ refers to the average of the layer with boundaries Z n 
and Z n+ i. If we write 
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then equations (A.l) and (A. 2) would become for I' energy points. 
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If we use " diamond scheme" [62] 
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and we write the equation (A.13) and (A.14) in the form of interaction prin- 
ciple. We obtain 
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where 



Xll X12 
X21 X22 



Yll Y12 
Y21 Y22 




+r[ s- ) (A. 19) 



Xll = M + -r(a + - 7]R ++ W ++ ^j 

X12 = --R + ~W + - 
2 

X21 = --tR- + W' + 
2 

X22 = M + \t(o- - R—W—^j 

Yll = M — -r(a + — — R ++ W ++ ) 
2 V 2 / 

Y12 = +-tR + -W + ~ 
2 

Y21 = +^tR- + W- + 

Y22 = M — \r[a~ - 

We can derive the transmission and reflection matrices and the source vectors, 
from the equation (A. 19) by the comparison with the interaction principle 
[62] and are given in the Appendix. 

To calculate the diffuse radiation ie., the specific intensities at each spacial 
point in the medium, 1^, we use the scheme described in [63] and [62]. We 
estimate 1^ at each point which are given by, 

I+ +1 = r(l,n+l)l; +1 + V+ (A.20) 
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I-=t(n,n+l)I- +1 + V~ H (A.21) 

where r(l, n+ 1), t(n, n+ 1) are the diffuse reflection and transmission matri- 
ces at the internal points. V + i , V~ x are the source vectors which represent 
the emission from the layer whose boundaries are Z n and Z n+1 together with 
the diffuse radiation from the rest of the medium. 

The quantity represents the reflected intensity in the direction of 
increasing r(r — > T, the total optical depth) while I~ represents the emergent 
intensity toward decreasing r — > 0. We can estimate these two quantities 
through the equations (A. 20) and (A.21) at any spacial point in the medium 
in the form of angular distribution. 

The transmission and reflection operators as are given by, 

t(n+l,n) = G + -[A+A+ + g + ~ + g~+] 

r(n + 1, n) = G~ + g- + [E + A+A+] 
t(n, n+l) = G-+[A-A- + g~ + g + 1 
r(n, n+l) = G +_ g +_ [E + A~A~] 
Where E is identity matrix and 

G + - = [E-g+-g- + ]- 1 
G- + = [E-g- + g+-]- 1 
g + - = \t±+ 
g -+ = - T A Y 

A+ = [M+ irZH 1 ; Z = a - 77RW; Y = r/RW 

and the source vectros are, 

S+ = rG + -[A + + g+-A-]S 

and 

E- = rG-+[A- + g-+A+]S 
S being the source function. 
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